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Abstract. A technique is proposed for constructing three-dimensional solutions 
comlplying with the self-consistent magnetohydrostatic (MHS) equations and with 
observations along the line of sight of the magnetic field at the photosphere. The 
^■y technique is a generalization of a paper (Rudenko, 1999) to a potential-field approx- 

imation. The solution of the problem under consideration leans a representation 
of the magnetic field in terms of a scalar function, with its subsequent harmonic 
expansion in terms of a spherical functional basis that satisfies specified boundary 
conditions. A numerical realization of the proposed method is expected to permit a 
real-time modelling of three-dimensional magnetic field, temperature, pressure and 
density distributions. 
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In a paper (Rudenko, 1999) it was suggested a method to model a three- 
dimensional distribution of the above-photospheric magnetic field, base 
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1. Introduction 



o 

i-^ ■ on solving the Laplace boundary- value problem that corresponds to ac- 



tual one-component magnetic measurements at the photospheric level. 



Boundary conditions on the lower sphere were specified in (Rudenko, 
1999) by a distribution of the projection of the magnetic field onto a 
certain fixed direction corresponding to the line of sight at the time of 
measurement. In this paper, the same boundary conditions are used in 
solving a boundary- value problem for a certain class of self-consistent 
magnetohydrostatic (MHS) equations admitting of analytic solutions. 
Developing self-consistent physical models of chromospheric and coro- 
nal plasmas is one of the major tasks of theoretical astrophysics. Using 
the MHS approximation is justified when it is intended to study the 
current state of the medium or relatively slowly occurring macroscopic 
plasma processes. Also, because of the non-linearity of the MHS equa- 
tions, the problem of seeking all possible solutions is extremely compli- 
cated. Several approaches to constructing self-consistent MHS models 
are known to date, which permit analytic solutions to be determined. 
These approaches are based on some not fundamentally differing pre- 
liminary assumptions about the way in which the solutions behave. 
An outline of the approaches and research results on this subject may 
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be found in ((Low, 1985), (Low, 1991) , (Low, 1992), (Low, 1993a), 
(Low, 1993b), (Bogdan and Low, 1986), and (Neukirch, 1995)). The 
class of the models studied thus far is relatively narrow to pretend to 
a generality of the description of the real states of magnetic plasma. 
Nevertheless, these models, when used, would be expected to give at 
least qualitative insight into the possible correlations between the fun- 
damental physical plasma characteristics. There is no question that the 
MHS modelling is much higher in information content than the simplest 
modelling in the potential-field approximation of the magnetic field in 
extensive use to date. 

An important factor as to the extrapolation of the magnetic field 
and other physical characteristics in terms of a particular physical 
model is the ability to derive physical solutions that satisfy boundary 
conditions corresponding to actual measurements. Only then can we 
expect that the result of an extrapolation will represent the facts, and 
only then can we select - on the basis of many comparisons - the 
most adequate physical model for use in the study. Three methods 
of specifying boundary conditions are customarily used to extrapolate 
magnetic fields in the potential-field approximation: from the radial 
component B r (Neumann's classical boundary- value problem); from 
the component B\ (projection along the direction corresponding to the 
orientation of the line of sight at the time when each point crosses the 
surface of the central meridian:, and from the component B& (projec- 
tion along the line of sight corresponding to the time of measurement 
(Rudenko, 1999)). Only the second and third methods of specifying 
boundary conditions correspond to actual magnetic field observations. 
More specifically, the second method gives only an averaged picture 
of the magnetic configuration (to accumulate boundary conditions re- 
quires a full series of daily observations over the course of a complete 
rotation). The third method appears to be the most correct as only its 
application can provide the best fit to a current state of the magnetic 
field. As in the case of the potential-field approximation, of important 
significance is the statement and solution of a boundary- value problem 
in the MHS approximation with boundary conditions of the third type 
(from the component Bd). In this paper such a setting of the problem 
is realized for MHS models as suggested in (Neukirch, 1995) where 
analytic solutions were obtained for MHS models satisfying boundary 
conditions of the 1st type (from the component B r ). In this paper 
we consider the same MHS approximations, yet for different boundary 
conditions. We restrict the discussion only to the analytic solution of 
the problem formulated. It can be used in constructing a numerical 
algorithm for calculating a current state of solar plasma, and this issue 
will be the subject of further investigation. The problems associated 
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with scarcity of data from the averted side of the Sun will not be 
discussed here as they were covered in considerable detail in a previous 
paper (Rudenko, 1999) for the potential-field case. 



2. Setting up the problem, and the method of solution. 



Here we take advantage of the physical models whose substantiation 
and physical aspects are reported in (Neukirch, 1995). We restrict 
ourselves only to a brief formulation of the problem in the form of equa- 
tions. From the logical standpoint, the problem is conveniently broken 
up into three parts: one is the system of basic equations; the other 
involves imposing additional requirements on the function of current 
density (strictly speaking, a concrete definition of the physical model) ; 
and the third includes specifying the boundary conditions. The system 
of MHS equations is 



jxB-Vp - pVV> = 0, (1) 

V x B = noi, (2) 

V • B = 0. (3) 

Here B, j, p, pand ip are magnetic induction, current density, pressure, 
density, and gravitational potential, respectively. 

As shown in (Low, 1991) and (Neukirch, 1995), the current density 
that allows pressure and density from (1) to be expressed explicitly in 
terms of the magnetic field has the form: 

j = aB + VF(VV>-B,V>) x Vip, (4) 

where F is an arbitrary function of its arguments. As is (Neukirch, 
1995), we consider a simplified model: a = const, and 



r 3 /riA/r\ 2 
r 5-r-B=£(r)- 

Here the gravitational potential is specified as 



CM r A / CM\' 

/,(0)^ r .B=£(r)— r-B, £( r ) =*(^_J . (5) 



GM _ , GM 
$ = — V ^ = r > ( 6 ) 



with M being the mass of the central body. For £(r), the same variants 
as in ((Neukirch, 1995), (40)) can be considered 
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{k/r 2 , case I, 

1/r 2 - l/(r + a) 2 , case II, (7) 

1/r 2 + fcexp(— 2r/L) — q, case III, 

where fc, q, and L are constants. 

To formulate the boundary conditions we introduced in the direction 
selected a fixed unit vector d, corresponding to the direction of the line 
of sight, and require the fulfillment of the conditions 

d-B lr=R = B d (6,4>), (8) 

B|_ - 0. (9) 

Here Bd(9,4>) is some specified distribution on a spherical surface of 
radius R. 

Formulas (2-9) are sufficient for a full determination of the magnetic 
field in the region of space r > R. The hydrostatic balance equation (1) 
describes pressure and density variations associated with the magnetic 
field. At this point we shall not give the formulas describing the specific 
relation of these quantities to the magnetic field; their derivation and 
particular representation may be found in ((Neukirch, 1995), (25)-(31)). 
In ((Neukirch, 1995), see (24)), an analytic expression was obtained 
for the magnetic field in terms of the coefficient of spherical harmonic 
expansion of the scalar product (r • B). It can be demonstrated that 
the magnetic field satisfying (2-5) may be represented in a more general 
(than in (Neukirch, 1995)) form in terms of a single scalar function. It 
is easy to verify by mere substitution that the magnetic field that is 
expressed as 

B = |-rA + |L + V^-v)r|x (10) 
in terms of a certain scalar function x satisfying the equation 

(A + eL 2 + a 2 ) X = (11) 

satisfies equations (2-5). In (10, 11), a = ^oa, £ = /io£, and L is the 
angular moment operator (see (Neukirch, 1995)) 

L = -r x V. (12) 

i 

Note that equation (11) for \ ls identical to the equation for (r • B) 
in (Neukirch, 1995). Hence the initial problem can be reformulated 
as a boundary- value problem for satisfying (11) and the boundary 
conditions (8), (9). 
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To solve the boundary- value problem formulated we introduce a 
spherical system of coordinates (r, 6, <p) connected with the unit vector 
d (as done in (Rudenko, 1999)). Let the axis z that correspond to the 
selected spherical system of coordinate, be directed along d, and let 
the axis x, for definiteness, be taken to lie in the plane produced by 
the axis z and the heliographic axis z'. Finally, let the function \ be 
represented as an expansion into a series in terms of spherical harmonics 
of the function 

ooVL I 

V(r,9,^ = Rj2 E <f/i W H (cos 9) e™t (13) 

1=1 m=—l 

where C™ are the desired complex expansion coefficients 

cJ m = cT, (14) 

pm _ pm J ^J2ttw™ , Pf 1 is the Legendre function (see (Abramowitz 
and Stegun, 1964)), and 

l 

w ™ = J [P{ n (u)] 2 du = {l + l/2)' 1 (l + m)\/{l-m)\. (15) 

-l 

The expression (13) should be regarded as an exact solution if the 
summation with respect to the index I is made ad infinitum, or as a 
finite approximation if the summation is limited to the value of L (the 
main index of expansion) . Substitution of (13) into equation (11) gives 
an equation which must be satisfied by the function fi(r): 

(£2 ~ l V + !) " + 5 ) r ^ r ) = °- ( 16 ) 

For the sake of simplicity, we did not introduced one further summation 
into (13) - with respect to the indices of two independent solutions of 
equation (16). For the three cases of the dependences £(r) in (7), the 
solitons of equation (16) are expressed in terms of different variants of 
Bessel functions (for the explicit expression see in ((Neukirch, 1995), 
(40)-(44)). Out of them we may always choose linear combinations 
satisfying the upper boundary conditions (9). Moreover, for the case 
of a = we can impose on some surface \r\ = R w > R the boundary 
conditions of a radial magnetic field. Such conditions are used very 
widely in the technique for extrapolating a potential magnetic field in 
modelling the conditions on the source surface. As will be apparent from 
the ultimate expansion for a full magnetic field vector, such conditions 
can be obtained by satisfying the equality 
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j-r (ai/, (1) (r) + a 2 // 2) (r)) = 0, Mr) = ai/, (1) (r) + a 2 jf >(r). 

When such conditions cannot be satisfied in principle; however, 

if a is small, the condition (17) might still be used. In this case the 
field will be approximately a radial one. In any case, here fi(r) is 
taken to mean the solutions of equation (16), satisfying some of the 
upper boundary conditions. Explicit expressions for them can always be 
obtained. Upon substituting the expansion (13) into the expression for 
the magnetic field (10), we arrive at the expansion of the full magnetic 
field vector 



B(r)=i?^ 1^1(1+1) + % + V ± ^A ft?!™ 1 (cos e)e im <t>. 

1=1 m =-l ^ r T ' 

(18) 

This expansion coincides exactly in its form with the expansion in 
((Neukirch, 1995), (24)). It can be used to calculate the magnetic field 
once its coefficients are determined. Note that by virtue of the identity 
(r • B) = L 2 x, the coefficients c[™ are those of the initial expansion 
(r • B) that is represented by the expression ((Neukirch, 1995), (17)). 

To solve the boundary- value problem formulated, an explicit form of 
the expansion of the component (d • B) is required. Using the represen- 
tation of the operators L and Vj_ in a spherical system of coordinates 

T--( — 66 9 \ \7 -¥L— e<f> 8 MCA 
i \ e<p d6 smOd^J' V± r M + rsm9d<f>' 1 ' 

we obtain 



ooVL 1 ( A i j "1 

d-B = «E E crrid + !)-«+ (l-« a )£i|;r|^N(«)e 

1=1 m=-l L J 

(20) 

Here and e$ are unit vectors of the spherical system of coordinates, 
and u = cos 6. To determine the coefficients of the initial expansion we 
take advantage of the projection scheme from ((Rudenko, 1999)). We 

multiply the equality (20) by the factor pj m ^ (cos 9) e~ im ^ cos sin 0d<pd6, 
and integrate it on a sphere |r| = R over angular variables: 
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j / B d (6, ip)Pj ml (cos 9) e^™^ cos 6 sm 6 d<pd6 = 



ooVL 



2f J [« 2 fc(fc + 1) - iai?m«] p] m! («) P, W («) (21) 

z=i -1 L 



-u^.^f^pM {u)\du 



Here = [^ r /fe( r )]| _ R It is also assumed that fk(R) = 1, which 
can always be done. Using the orthogonality property of Legendre 
polynomials and the known recurrence relations for expressions of the 
form uP™ and (u 2 — \)dP™ /du ((Abramowitz and Stegun, 1964)), and 
introducing the designation for the left-hand side of the equality (21) 
bf 1 (these quantities will be referred to as the weight coefficients)we can 
bring the expression (21) into the form: 

h? = a m _ic?_ 2 + a^cf-i + <^T + <(<T +1 + a#c£ 2> (22) 
where 



m.l 



™T-l (Z-m-l)(Z+m)(Z-2) ,j , v 
y/w™ 2 w™ (2Z-3)(2i+l) V 1 91-2), 



m,Z _ /n /,',„ /_!!_?_ '■' 



a 



a +1 - tanm^i w ^ (2i+3) , 

m,Z _ (Z-m+l)(Z+m+2)(Z+3) /, „ s 

fl +2 ~ y/wfawF (2Z+5)(2Z+1) ^ + Z + 5Z+2J- 

(23) 

In view of the property of (14), when determining the coefficients cf 1 , it 
is sufficient to use equations (22) for positive m only. The equalities of 
(22) can be broken up into m independent systems of equations, each 
of which can be represented in the matrix form (for a finite expansion 
to L) 

A™C? = B™, (24) 

where 
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,m,2j+m~2 _|_ fii 
m,2j+m—2 



•i m,2j+m-2 
] j a 




,rn 



■2j+m-2i D i 



Br = 



'2i+m-2i 1 - 



i,j < {L-m)/2 + l,6 i j 



M = 3 
(25) 



Hence, by determining the coefficients c from equation (24), we solve 
the problem formulated. For each value of m, the matrices A m in (24) 
have a pentadiagonal form. Since the selected values of a are always 
small and the behavior of the elements on even diagonals is similar 
to the behavior of the respective elements of the matrix derived in 
((Rudenko, 1999)), one might expect that calculations of this vector 
equation will also present no problems. 



The formulas obtained in this study for determining the expansion 
coefficients of the magnetostatic solution can be used in further inves- 
tigation in calculations using real daily magnetograms. As shown in 
((Rudenko, 1999)), this technique can only work well with sufficiently 
high resolution data, such as magnetograms from Kitt Peak Observa- 
tory. Technically, the realization of the numerical code of the solution 
outlined above is not too different in its complexity from the code of 
calculating an analogous problem for the potential-field model. One 
might expect good results from extrapolating solar plasma parameters 
throughout the visible disk at the time of observation, to the heights 
of about the solar radius ((Rudenko, 1999)). 
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